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The hole spectral function of the t-J model on a two-dimensional semi-infinite lattice is calculated 
using the spin- wave and noncrossing approximations. In the case of small hole concentration and 
strong correlations, t ^ J, several near-boundary site rows appear to be depleted of holes. The 
reason for this depletion is a deformation of the magnon cloud, which surrounds the hole, near the 
boundary. The hole depletion in the boundary region leads to a more complicated spectral function 
in the boundary row in comparison with its bulk shape. 
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I. INTRODUCTION 

In recent years, an active interest is taken in the 
electronic properties of heterostructures and surfaces of 
strongly correlated materials.— Looking for new effects 
and their possible applications a wide variety of sys- 
tems has been investigated both experimentally and the- 
oretically. Theoretical studies of charge excitations near 
the crystal boundary have been carried out mainly in 
the framework of the two- (2D) and three-dimensional 
(3D) Hubbard model. For this purpose different ap- 
proximate methods have been used, including the slave 
boson method,^ the perturbation theory^ and the dy- 
namical mean-field theory.i^ In these works, the case of 
half-filling was considered, when strong electron correla- 
tions cause the antiferromagnetic ordering of the crystal.^ 
However, approximations used in the mentioned works 
did not take into account the ordering and the interaction 
of electrons with respective magnetic excitations. One of 
the results obtained in Refs. y-la for uniform model pa- 
rameters is that on the surface layer the quasiparticle 
weight is smaller than the bulk value. The reason is a re- 
duced surface coordination number which implies a lower 
kinetic energy and consequently effectively stronger cor- 
relation effects at the surface. 

It is knownl that in the case of strong electron cor- 
relations the interaction with magnetic excitations plays 
an important role in the formation of the low-frequency 
dispersion of charge carriers. Therefore, peculiarities of 
these excitations in the near-boundary region may have a 
significant impact upon the properties of electrons here. 
The magnetic excitations are described by the quantum 
Heisenberg modeli^ The influence of boundaries on its 
spectrum and observables has been studied in two^""— 
and three^^ dimensions. In particular it was shown that 
absolute values of the nearest-neighbor spin correlations 
near the boundary exceed the bulk value. In Refs. [lllfT^ 
this result was related to the peculiar spectrum of the 
semi-infinite d-dimensional antiferromagnet. The spec- 
trum involves d-dimensional bulk modes ~ standing spin 
waves - and a (d— l)-dimensional mode of boundary spin 
waves. These latter excitations eject the bulk excitations 
from the near-boundary region. Thus the antiferromag- 
net appears to be divided into two regions with different 



dominant spin excitations. Charge carriers in the near- 
boundary region and deep within the crystal appear to be 
in different spin-excitation environments, that inevitably 
leads to a dissimilarity in properties of these carriers. An- 
other effect which can contribute to this difference is a 
smaller number of spin bonds destroyed by charge car- 
riers near the boundary in comparison with the bulk. 
As will be seen below, this leads to an attraction of the 
quasiparticles to the boundary. 



To answer the question on how the above-mentioned 
factors influence the distribution of charge carriers near 
the boundary we consider the 2D t-J model, having in 
mind a Cu-0 plane of hole-doped cuprate perovskites. To 
calculate the hole Green's function in the case of strong 
correlation, t ^ J {t and J arc the nearest neighbor 
hopping and exchange constants) we apply the spin-wave 
and noncrossing approximations which were successfully 
used for unbounded crystals J^ii^ An obtained self-energy 
equation for the hole Green's function is solved by itera- 
tions. Since the translation invariance is violated in the 
X direction perpendicular to the boundary, the spectral 
function, apart from the frequency and the y component 
of the wave vector, depends on the x coordinates of site 
rows for which the function is considered. With a change 
of the X coordinates from (the boundary) deep into the 
crystal the intensity is redistributed in the function and 
the main maximum is enhanced and shifted to lower fre- 
quencies. Such behavior of the maximum indicates that 
the near-boundary region is depleted of holes at low hole 
concentrations. The appearance of this depletion is con- 
nected with the character of hole excitations. They are 
spin polarons in which a hole is surrounded by a cloud 
of magnons. Near the boundary, this cloud is deformed, 
which leads to an energy loss and to the observed shift of 
the main maximum to higher energies. Another conse- 
quence of the frequency-separated maxima in neighbor- 
ing rows is a more complicated structure of the boundary 
spectral function in comparison with its bulk counter- 
part. The reason is a replica of a stronger maximum 
from the second row which is also seen in the boundary 
spectral function. 
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II. MAIN FORMULAS 

Our starting point is the Hubbard Hamiltonian on a 
square semi-infinite lattice. We consider an idealized 
boundary which is located along the y crystallographic 
axis. The variation of the lattice spacing and model pa- 
rameters near the boundary is neglected. The Hamilto- 
nian reads 



lyCr l^>0 



where a; i^^- is the electron annihilation operator, ly and 
Lj: are the vector components labeling sites of the crystal, 
which is located at > 0, cr = ±1 is the spin projection, 
6 ~ ±1, the lattice spacing is set as the unit of length, 
U is the Hubbard on-site repulsion, ni i „ = a] , „ai i a, 

i I y X lyi^u y X 

and fi is the chemical potential. Only the hopping be- 
tween nearest neighbor sites t is taken into account in 
Eq. (P. 

In the case of strong electron correlations, U :s> t, and 
an electron filling less than half-filling Hamiltonian ([T]) 
can be reduced to the Hamiltonian of the t-J model using 
the known unitary transformation^^ Htj = e^Hffe^^ 
with 



u 



ly6(T Ix >0 

t 
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^lyl^ ^lyl:. + i ^lyl^^ly,lr. + l I ' 



where the Hubbard operator: 
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aiyl^a = + ^^^„^J^ (^lyl^a = + '^^fylf 

were introduced. Up to the terms of the second order in 
the transformed Hamiltonian reads 

ly6(T Ij^ >0 

+t V V fxfO; +Xfi +1 

+ J X! \^ly + l,l^^lyl:. +^ly,l.v + l^ 



5/ / 



+ - ■^^ly,l^ + l'^lyU 



(2) 



where S;^/_^ is the spin-i operator, J = and ni^i^ = 
X](T "ia'xCT = 1 — -^lyi^ the considered approximation 
in which terms containing doubly occupied site states 
are neglected. In Eq. ([2]), we neglected also terms pro- 
portional to J which describe an assistant hole hopping 
(three-site terms), as it is frequently done in the consid- 
eration of the t-J model. 

Further simplifications of the model Hamiltonian are 
connected with the spin-wave approximation which in 
application to the t-J model was shown to give results 
in good agreement with exact diagonalization.-'^'^'^^ In 
the case of low doping and zero temperature the crys- 
tal has the long-range antiferromagnetic ordering and 
the simplest version of the spin-wave approximation 
can be applied using the following Holstein-Primakoff 
representationiii 

S^yl^ - Plyl^VlylMyl. + P^ylAyl.^lyl^^^ 

where the spin-wave operators 6/^;^ and b\ j satisfy the 
Boson commutation relations and 



In the considered antiferromagnetic background the hole 
creation operator is defined as 



=y 

lyix / J I'y I'X t y I'X 



Using this definition and Eq. ([3)) in Hamiltonian ([2]) and 
leaving terms up to the second order in the spin-wave 
operators we get 



H = 



tJ2 J2 K+s,iMyi^ (biyi^ + bi^s,i^) 



+t 



J2 [K^l^+i'^lyl^{\l^+blj^ + l 

ly,l^>0 

+ KiAyX^+l {K^.+l+^\yl^) 
+ HAF ~ ^ E {'^ly+l,lx'^lylX + l^lyAX+ll^lylx) 



where i>i 



= h\yi.Ki^ and 



(4) 



haf ^ J \^(^-ls^^^)bll^yl. 

ly.la:>0 L ^ ^ 



(5) 
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is the Hamiltonian of the 2D semi-infinite Heisenberg 
antiferromagnet in the spin-wave approximation. In 
Eq. (jll) , some constant terms were omitted and the term 
I J was added to the chemical potential. 

The next to last term in the right-hand side of Eq. (|4]) 
describes an attraction of a hole to the boundary. It 
originates from terms of Hamiltonian ([2]) which contain 
z components of spins and occupation numbers on neigh- 
boring sites. In the antiferromagnetic state, these terms 
give the energy gain equal to ^ for each nearest-neighbor 
bond. In the 2D case a hole destroys 4 such bonds deep 
inside the crystal and 3 bonds on the boundary. Thus, 
for a hole it is energetically more favorable to reside at 
the boundary. 

Refusing the constraint Za: > and carrying out 
the Fourier transformation over the space coordinates, 
Eq. (|4]) is reduced to the spin-wave Hamiltonian on an 
unbounded lattice, used in Refs. [Tslfl^ and in a lot of 
subsequent works. 

Considering the case of a low hole doping, in Hamilto- 
nian we shall neglect terms containing two hole oc- 
cupation operators on neighboring sites. Our aim is to 
calculate the hole Green's function 



where the angular brackets denote the statistical aver- 
aging, T is the time-ordering operator that arranges 
other operators from right to left in ascending order 
of times r, /ife ; is the Fourier transform oi hi i , and 



^kyi^{T) — hkyi^e^^ ■ For this calculation, we use 
the diagram technique with the expansion in powers of 
t, also in full analogy with what was done for the un- 
bounded crystalJ^i^ For this latter case, the self-energy 
equation was obtained in the noncrossing (Born) approxi- 
mation in which diagrams with intersecting magnon lines 
were neglected. It was shown that results obtained in this 
approximation are in good agreement with data of exact 
diagonalization. Therefore, we also use this approxima- 
tion and obtain the following self-energy equation: 



■xY.{kynlxilx2)G{kynlx2l'x), 

kyi/ s,s' — — 1 

xG{ky — k'y,n ~ I'Jx + s, l'^ + s') 
X gky~k'^,s9kys'Di2{k'yi^lx,l'^ + s') 
+9ky-k'^ .sQky -k'^,s' Dii{k'yVlJ'^ ) 

+ 9kySgkyS'D22{k'yiy, 

+9kysgky-k'^,s'D2i{k'yV, Ix + s,Q 



(6) 



(7) 



where n and arc shorthand symbols for the Matsubara 
frequencies a;„ = (2n -I- 1)ttT and 2vnT ^ respec- 

tively, T is the temperature, 

ei^ = fi — ^Si^o with the last term taking into account 
the attraction of a hole to the boundary, N is the number 
of sites in the y direction. 



2t coa{ky), s — 0, 
t, s = ±1, 



and Dijlkyvlxl'x) is the Fourier transforms of the compo- 
nents of the matrix magnon Green's function 



D{kyTlX) 

Bk^i^ 



rSkyMBlu 

bk^i^ 



(8) 



Equation ([7]) describes the contribution of the sunrise 
diagram. In the case of the semi-infinite crystal there is 
also a nonzero contribution of the bubble diagram, which 
vanishes in an unbounded crystal. One can show, how- 
ever, that in the semi-infinite crystal in the case of small 
hole concentrations the contribution of the bubble dia- 
gram is also negligibly small. Indeed, this term contains 
the multiplier 



^9kys{hk 



(9) 



where the mean value can be expressed through the re- 
tarded hole Green's function G{kyUjlxl'x) as 



duj lmG{kyUjlxl'x) 



1 



.IT 



For T = the integration is carried out over unoccupied 
states. For small hole concentrations these are in fact all 
states and therefore one can rewrite the above integral 
as 



/oo 
— lmG{kyU}lxlx) 
-oo 



Thus, for s = ±1 the multiplier © is small because 
the mean value is negligible, while for s = it is small 
because g^yO ex. cos(fcy) and the sum over ky is negligible. 

Let us switch from the Matsubara Green's functions to 
the real-frequency retarded Green's functions. It can be 
done using the following relation between these functions 
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doj Im [Dij{kyU}lJ'x) + Dji{kyLol'Jx)] - »Re [D,j{kyOjlJ'x) - DjjjkyUjl'Jx)] 
27r uj — iu),. 



(10) 



The relation can be verified using the spectral representations. An analogous relation can be written for the hole 
Green's functions. From Eqs. ([6]), ([7]) and equations given below one can see that 

Dij{kyUjlJ'x) = Dji{kyUjl'Jx), G{kyLolJ'x) = G{kyUjl'Jx). 



Thus, only imaginary parts of the retarded Green's functions appear in Eq. (|T0|) . Substituting these representations 
into self-energy and carrying out the summation over v we find 



1 ^ /"OO T / 

lmJ:ikyU;lX) = E 0{L + s)9{l', + s') J^l^^ky ~ k'y., 

k' s,s' = -l ^ 



UJ — Oj', Ix + S,l'x+ s) 



9ky-k'^,sgkys'lmDi2{k'yUj'lx,l'x + s') + gky-k'^,s9ky-k'^,s''imDii{k'yUj'lxQ 
+9kys9kyA^D22{k'yUj' ,lx + s,l'^ + s') + gkys9ky-k'y,A^D2i{k'^uj' ,lx + s,Q 



(11) 



with nj^(w) = (e'^/^ + l) ^ and ns(w) = (e"/^ - l) \ 
The real part of self-energy PT|) can be calculated from 
the Kramers-Kronig relation. Self-energy equation (jG)) is 
transformed to real frequencies by the substitution za;„ — > 
cj + i?7, 77 — > +0. 

In the considered case of small hole concentrations we 
can neglect the influence of holes on magnon Green's 
function ^ and use its value for the undoped case de- 
scribed by Hamiltonian ([5]). In this case Green's function 
read o^-^-^^ 

b{kyUjlj'x) = D^"\kyUX) - ^D^'^HkyUjlxO) 



-D^°\kyuj00) 



where / is a 2 x 2 identity matrix, 



dkx sm[kxilx 



UJ — i?k + ir] UJ + _Ek + 

2 



D'^°\kyuj0l'x), 
(12) 

l)]sm[kx{l'x + l)] 
(13) 



Q , Uk«k 



k = {kx,ky), i?k 2J-\/l — 7^ is the bulk magnon en- 
ergy, 7k = 5 [cos(fc3,) -f cos(fcj^)], and 



Uk 



1 
2 
1 

fk = — 

2 



1 " 7k 

l + 7k 



1 + 7k 



1 


+ 7k 




-7k 


1 


+ 7k 




-7k 



In Eq. (IT^ . the first term in the right-hand side de- 
scribes the bulk modes - the standing spin waves, while 
the second term is connected with the boundary spin 
waves. Their peak dominates in the spectral intensity 
—lv[\D{kyUjlxlx) for /a; = 0, 1 and practically disappears 
in site rows more distant from the boundaryj ^^'^^ 

It is instructive to elucidate how the equations ob- 
tained above are transformed to the form for an un- 
bounded crystal with distance from the boundary. In 
Eq. (fT2|) . the second term in the right-hand side becomes 
negligibly small if at least one of the coordinates Ix or 
I'x is larger than 2. Green's function D'-^^kyUjlxl'x), to 
which D{kyUilxl'x) is reduced for such x coordinates, con- 
tains the multiplier sm.\kx{lx + 1)] sin[fca;(Z^ + 1)] in its 
integrand [see Eq. (|13p ]. If in this multiplier the sines 
are replaced by their representation through exponential 
functions, one can realize that terms with the same signs 
of exponents are small for large Ix or since the respec- 
tive exponential functions rapidly oscillate. Remaining 
terms depend only on the difference Ix — I'x as it must for 
the unbounded crystal. It can be shown that these terms 
are identical to Green's function for this case. Since the 
magnon Green's function defines the form of the hole 
Green's function, one can expect that the latter also be- 
comes close to its unbounded form when at least one 
of the X coordinates is large. Taking this into account, 
from Eqs. (O and pT|) - (fT5|) after the Fourier transfor- 
mation we obtain equations for the unbounded crystal of 
Refs. [illll 

The above discussion allows us to transform Eq. ^ 
into a more tractable form. Let us rewrite it as 



Zl^)5l^in-ll{kyUjlxl'l)\ GikyUjl'X) 
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= 5i^i,^ + M{kywlX), (14) 

where 

In Eq. ([14]), we assume that the coordinates Ix and l'^ are 
restricted within the range [0,^2;^]- The parameter l^m 
is expected to be large enough for substituting the self- 
energy and Green's function in Eq. (jlSp by their values 
in an unbounded crystal, in compliance with the above 
discussion. At the same time Ixm can be chosen to be 
small enough for the inversion of the matrix in the left- 
hand side of Eq. (fT4| would not lead to time-consuming 
calculations. 



III. RESULTS AND DISCUSSION 

In the below calculations we set T = and J/t = 0.2. 
The latter ratio of parameters corresponds to hole-doped 
cuprates .^^1^^ Equations (IT4t and (|15p were solved by it- 
erations for Ixm = 4, using as the starting function for 
G{kyUjlxl'x) Green's function of an unbounded crystal. To 
ensure the convergence of the iteration procedure an ar- 
tificial broadening was introduced by substituting uj with 
u! + ir], rj = O.Obt, in Eq. (fM)) . The chemical potential 
fi was chosen so that the frequency a; = 0, which sepa- 
rates occupied and unoccupied states, was located in the 
low-frequency tail of the spectral function 

A{kyLolx) ~ —lmG{kyUjlxlx)- 

This ensures a low hole concentration expected in the 
derivation of the above formulas. ■^'^ 

The spectral function gives the density of states pro- 
jected on states of the row Ix- A typical example of 
this quantity obtained in the course of the calculations is 
shown in Fig. [1] Besides the spectral function for near- 
boundary rows. Fig. [1] contains also the spectral function 
of the unbounded crystal, 

1 r 

Ab{kyUj, Ix - I'x) = -TT- / dkx cos [kx {Ix - I'x)] 

xliJiGbikykxid), 

which is given for comparison. Due to the translation 
symmetry this function depends only on the difference 
Ix — I'x and for the considered case Ix — I'x its last argu- 
ment is zero. In shape this function resembles spectral 
functions obtained for a fixed wave vector kx in an un- 
bounded crystalf^^^ However, the maxima in Fig. [1] are 
somewhat broadened in comparison with these functions 
due to the integration over kx in the above formula. As 
would be expected, the spectrum in the boundary row 
Ix = differs most greatly from Ai,{kyUjO). From the fig- 
ure one can see how the spectrum is transformed, gradu- 
ally approaching to the spectrum of an unbounded crys- 
tal, with distance from the boundary. In the scale of 
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FIG. 1: (Color online) The spectral function A(kyLjlx) for 
ky = 0, Ix — 0, 1, 4 and in an unbounded crystal. 
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FIG. 2: (Color online) The spectral function Aikyiolx) in the 
vicinity of the main maximum for fcy = (a), 7r/4 (b) and 
7r/2 (c) in the rows i^, = 0, 1, 4 and in an unbounded crystal. 



Fig. [T] already the spectrum in the fifth row (Ix = 4) is 
barely distinguishable from Af,(fcywO). 

The vicinity of the main maximum of the spectral func- 
tion is shown in Fig. [2] for several wave vectors. From 
formulas of the previous section it can be shown that 
Ai^iyiolx) = A(tx — ky,ujlx). Besides, in accord with the 
symmetry A{kyUjlx) — A(~ky,ujlx). From these figure 



6 



and equations one can see that the spectral maxima are 
shifted to higher frequencies on approaching the bound- 
ary for all wave vectors ky. In accordance with this the 
low-frequency tails of the maxima become weaker with 
decreasing l^- Since the concentration of holes in a row 
x{lx) is determined by this tail, 

ky 

u J-oo ~ 

one can conclude that the concentration decreases mono- 
tonically on approaching the boundary - near-boundary 
rows are depleted of holes. 

To elucidate a formation mechanism of this hole de- 
pletion layer let us first consider the role of two above- 
mentioned factors, which can influence the population 
of holes in the near-boundary region - the attraction of 
a hole to the boundary and the near-boundary magnon 
mode. For the chosen chemical potential, the energy of 
an immobile hole is equal to 6t. As seen in Figs. [T] and [21 
by virtue of the interactions there is the energy gain equal 
approximately to 2.5i in states corresponding to the main 
maximum. Therefore, the attraction which is of the or- 
der of J ^ t plays practically no role in the hole distri- 
bution. This conclusion is confirmed by calculations - 
omitting the attraction term is barely perceptible in the 
shape and location of maxima of the spectral function. 
The contribution of the near-boundary magnon mode can 
be evaluated from Fig. [31 The spectral function without 
this mode was calculated with the magnon Green's func- 
tion instead of the full function As seen from 
the figure, the near-boundary mode makes its contribu- 
tion in the location of the maximum and the intensity 
redistribution. This contribution is especially detectable 
for the boundary row. However, with this mode or with- 
out it the main maxima in the near-boundary rows have 
higher frequencies than deep inside the crystal. Conse- 
quently, the near-boundary mode does not play the main 
role in the formation of the hole depletion layer. 

To understand the appearance of the depletion layer 
let us remind that in the considered model holes are spin 
polaronsi ^^i^^ Due to the antiferromagnetic background 
a hole can move over the lattice only with the emission 
and absorption of magnons, as it is seen from Hamilto- 
nian dJ]). As a consequence the hole is surrounded by a 
cloud of magnons. Without spins the maximum energy 
gain which a moving hole can achieve in comparison with 
an immobile quasiparticle is At - the difference between 
the lowest energy in the 2D nearest-neighbor band and 
its center of mass. In the antiferromagnetic lattice this 
gain is decreased by the energy consumption for the dis- 
tortion of the magnetic order around the hole. For the 
ratio J/t = 0.2 the energy gain is reduced approximately 
to 2.5t (see the above figures). This energy gain is still 
comparable with the maximal possible value At. Notice 
that at the same time the spin polaron bandwidth is of 
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FIG. 3: (Color online) The spectral function A{kyUjlx) in the 
vicinity of the main maximum with taking into account the 
boundary magnon mode (black solid lines) and without it 
(red dashed lines) for ky = 7r/4, Ix = (a) and 1 (b). Blue 
dash-dotted lines corresponds to At{kyijO). 



the order of J for low doping, which is much smaller in 
comparison with the energy gain.— This large energy 
gain complicates the formation of ferrous - ferromagnet- 
ically ordered regions around holes - and stripes in the 
t-J model. Only for very small ratios J/t the gain in 
the hole kinetic energy in the ferromagnetic region be- 
comes large enough to stabilize ferrous i^""— Away from 
the boundary the magnon cloud has the symmetry de- 
termined by the group of the hole wave vector and this 
symmetry ensures the lowest energy of the spin polaron. 
Near the boundary, the cloud is distorted, which lowers 
the symmetry and inevitably leads to a growth of the 
energy. It is the mechanism of the depletion layer forma- 
tion in the considered model. The depth of the row in 
which the location of the main maximum coincides with 
that in an unbounded crystal gives an estimate of the 
magnon cloud size. In our case, its radius is equal to 4 
lattice spacings. 

Notice that as in Refs. in our case the decrease 
in the spectral intensity of the main maximum in the 
boundary row is connected with the reduced boundary 
coordination number. In the mentioned works this leads 
to an effective strengthening of the Hubbard repulsion on 
the boundary, while in our case to the deformation of the 
magnon cloud around the hole in the spin polaron. 

Closer inspection of the obtained spectral functions 
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FIG. 4: (Color online) The spectral function A{kyLjlx) in the 
vicinity of the main maximum for Ix ~ and ky ranging from 
(the bottom curve) to 7r/2 (the upper curve) with the step 
7r/20. For better visibility curves with larger ky are shifted 
upward with respect to curves with smaller wave vectors. 



FIG. 5: (Color online) The same as Fig.[4l but for the spectral 
function of an unbounded crystal Ab{kyUjO). 



shows that the main maximum for the boundary row 
has a more complicated structure than spectra for 1^.-2 
to 4 and for the unbounded crystal. Figures |4] and [5] 
demonstrate this difference. The evolution of maxima in 
rows = 2 to 4 are similar to that shown in Fig. [SJ 
while for — I the spectrum has some features of the 
boundary row. This result demonstrates how the deeper 
near-boundary regions approach in their properties to the 
unbounded crystal, which states are characterized by a 
2D wave vector. As known, ^-li^ the energy of these states 
has a minimum at the points (±^, ±-|) and their disper- 
sion is weak on the boundary of the magnetic Brillouin 
zone, which is composed of segments (0,±7r) — (±7r, 0). 
The states near these segments make the main contribu- 
tion into the maximum in Fig. [5]- for ky = wave vectors 
of these states lie near (±7r,0), while for ky — ^ these 
wave vectors are from the neighborhood of (±f , f) ■ The 
change in the location of the maximum when ky grows 
from to |- in Fig. [5] reflects the mentioned weak dis- 
persion of the states along the boundary of the magnetic 
Brillouin zone. The shoulder, which approaches the max- 
imum from high frequencies, is mainly connected with 



states from the vicinity of the axes and the boundary of 
the Brillouin zone - for ky — ^ these states have wave 
vectors near (0, ^) and (±7r, A similar high-frequency 
shoulder is observed in Fig. SI However, in contrast to 
the deeper rows, the main maximum in the boundary row 
has also a low-frequency shoulder which is best seen for 
small ky. As follows from Fig. [2lja), the locations of this 
latter shoulder is close to the position of the maximum 
in the row — 1. Indeed, in the considered system two 
neighboring rows have maxima, which are shifted in the 
frequency scale relative to each other. Since in accord 
with the formulas of the previous section the spectral 
function of a given row is connected with functions in 
neighboring rows, one can expect that a replica of the 
more intensive maximum for Ix = 1 will be seen in the 
boundary row. In the present case this replica looks like 
the low-frequency shoulder of the main maximum. Thus, 
a more complicated character of the boundary spectra is 
connected with the replica of the maximum of the under- 
lying row. Notice that this replica is an attendant effect 
of the hole depletion in the near-boundary region. 
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IV. CONCLUSION 

Our calculations referred to a 2D crystal. From the 
similarity of the 2D and 3D magnon spectr a^^'^^ we can 
expect to obtain analogous results for charge carriers in a 
3D crystal with strong electron correlations. From these 
results, the conclusion can be drawn that the surface 
electronic structure, which is tested by the photoelec- 
tron and tunnel spectroscopies, even in the considered 
case of the idealized surface may essentially differ from 
the bulk spectrum. The discrepancies between the pho- 
toemission data of a number of transition-metal oxides 
and calculated bulk spectra were interpreted similarly in 
Refs. 24.25. 

Comparing results obtained in the semi-infinite Hub- 
bard^r^ and t-J models, we find some similar features. In 
spite of the differences of models and computation meth- 
ods, in both models for uniform parameters the quasi- 
particle weight is lowered, while the intensity of the high- 
energy part of the spectrum grows^ at the boundary. The 
reason for this intensity redistribution is similar - it is 
a reduced coordination number at the boundary, which 
leads to an effective strengthening of the on-site repul- 
sion in the Hubbard model and to the deformation of the 
magnon cloud around the hole in the spin polaron in the 
t-J model. 



In summary, in the present article, we investigated the 
spectral function of the 2D t-J model on a semi-infinite 
lattice. The limit of strong electron correlations, t ^ J, 
and the case of low hole concentrations were considered. 
For this investigation, we used the spin-wave approxi- 
mation and the diagram technique with the non-crossing 
approximation. The obtained self-energy equations were 
solved by iterations, and we could trace the variation of 
the spectral function with distance from the boundary. 
Already in the fifth row the spectral function nearly co- 
incided with its bulk counterpart. It was shown that the 
near-boundary region of the crystal is depleted of holes. 
The reason is the deformation of a magnon cloud around 
a hole in this region, which is accompanied by energy 
losses. The hole depletion is reflected in dissimilar lo- 
cations and intensities of the main spectral maxima for 
different site rows near the boundary. As a consequence 
a replica of a maximum in the second row is seen in the 
boundary spectral function. This results in its more com- 
plicated shape in comparison with the bulk spectrum. 
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